arXiv:l502.00606vl [math.DS] 2 Feb 2015 


NO HYPERBOLIC PANTS FOR THE 4-BODY 

PROBLEM 

CONNOR JACKMAN AND RICHARD MONTGOMERY 


Abstract. The iV-body problem with a 1/r 2 potential has, in 
addition to translation and rotational symmetry, an effective scale 
symmetry which allows its zero energy flow to be reduced to a 
geodesic flow on complex projective N — 2-space, minus a hyper¬ 
plane arrangement. When N = 3 we get a geodesic flow on the 
two-sphere minus three points. If, in addition we assume that the 
three masses are equal, then it was proved in (IJ that the corre¬ 
sponding metric is hyperbolic: its Gaussian curvature is negative 
except at two points. Does the negative curvature property persist 
for N = 4, that is, in the equal mass 1/r 2 4-body problem? Here 
we prove ‘no’ by computing that the corresponding Riemannian 
metric in this N = 4 case has positive sectional curvature at some 
two-planes. This ‘no’ answer dashes hopes of naively extending 
hyperbolicity from N = 3 to N > 3. 


1. Introduction 

In [I] it was shown that the reduced Jacobi-Maupertuis metric for a 
certain three-body problem had negative Gaussian curvature (except at 
two points where it is zero). This hyperbolicity led to deep dynamical 
consequences. Does hyperbolicity, i.e. curvature negativity, persist 
for the analogous iV-body problem, N > 3? No. We show that the 
analogous reduced 4-body problem with its metric has two-planes at 
which the sectional curvature is positive. 

The iV-body problem in question has equal masses and the inverse 
cube law attractive force between bodies. 

2. Set-up 

Identify the complex number line C with the Euclidean plane R 2 .Then 
the planar iV-body problem has configuration space C N \ A. Here A is 
the “fat diagonal” consisting of all collisions: A = {q = (gq, g 2 , ■ ■ •, Qn) £ 
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CA : q t = qj for some pair i 7 j}. The quotient of CA \ A by transla¬ 
tions and rotations is the “reduced iV-body configuration space”: 

C N = Y n x M+ Y n = CP^" 2 \ PA 

where CP W ~ 2 is the projectivization of the center of mass subspace 
CA -1 = {q e CA : Yriipp = 0} and PA C CP^ 2 is the projectivization 
of A fl CA -1 . The M + factor records the overall scale of the planar N- 
gon and is coordinatized by \fl with 7 = Emj|gj| 2 being the total 
moment of inertia about the center of mass. Y N is the moduli space 
of oriented similarity classes of non-collision TV-gons and will be called 
“shape space.” 

The following considerations reduce the zero angular momentum, 
zero energy iV-body problem to a geodesic flow on shape space Y/v, 
provided the potential V is homogeneous of degree —2. If A is ho¬ 
mogeneous of degree —a then the virial identity, also known as the 
Lagrange-Jacobi identity, asserts that along solutions of energy 77 we 
have 7 = 477 — (4 — 2a)V which implies that the only case in which 
we can generally guarantee that I = 0 is when a = 2 and 77 = 0. If 
in addition 7 = 0 then solutions lie on constant levels of I. Now we 
recall the Jacobi-Maupertuis [JM] reformulation of mechanics which 
asserts that the solutions to Newton’s equations at energy 77 are, af¬ 
ter a time reparameterization, precisely the geodesic equations for the 
Jacobi-Maupertuis metric 

ds 2 JM = 2(77 - V)ds 2 

on the Hill region {77 — V > 0} C CA \ A with ds 2 the mass metric. 
We are interested in the case 77 = 0, —V > 0 with V homogeneous of 
degree —2, in which case the Hill region is all of CA \ A and 

ds 2 JM = Uds 2 , U = -A 

The case of prime interest to us is 
(1) U = — A = N h//,///,7/' 2 

where the sum is over all distinct pairs ij. This U, and hence the JM 
metric, is invariant under rotations and translations. Qnotienting first 
by translations we take representatives in the totally geodesic center of 
mass zero subspace CA -1 , which reduces the dynamics to geodesics of 
the metric ds 2 JM \cN~i on CA -1 . Moreover, ds 2 JM \cN-i is also invariant 
under scaling since the homogeneities of U and the Euclidean mass 
metric ds 2 on CA -1 cancel. Thus the JM metric admits the group 
G = C* of rotations and scalings as an isometry group. Now Y N is 
the quotient space: Y N = (CA -1 \ A )/G = CP jV ~ 2 \A. (By abuse 
of notation, we continue to use the symbol A to denote the image of 
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the collision locus A under projectivization and intersection.) Insisting 
that the quotient map n : \ A —>• Y/v is a Riemannian submersion 

induces a metric on Y/v- Recall that this means that we can define the 
metric on Y n by isometrically identifying the tangent space to Y n at 
a point p with the orthogonal complement (relative to ds 2 JM or ds 2 , 
and at any point lying over p in C^ -1 ) to the G-orbit that corresponds 
to that point. These orthogonality conditions are equivalent to the 
conditions that the linear momentum, angular momentum, and ‘scale 
momentum’ / are all zero. To summarize, by using the JM metric 
and forming the Riemannian quotient, the zero-angular momenetum, 
zero energy 1/r 2 IV-body problem becomes equivalent to the problem 
of hireling geodesics for the metric defined by Riemannian submersion 
on Y n . 

Remark. The metric quotient procedure just described realizes the 
Marsden-Weinstein symplectic reduced space of T*(C N \ A) by the 
action of translations, rotations and scalings, CxiC*, at momenta values 
0, together with the IV-body reduced Hamiltonian how, but valid only 
at zero energy. 

Remark This metric on Yn can be expressed as Uds 2 FS where ds 2 FS 
is the usual Fubini-Study metric on CP jV ~ 2 . 

Remark For the standard 1/r 2 potential of (eq. [T|) this metric on 
Yv is complete, infinite volume. 

The collincar IV-body problem defines a totally geodesic submanifold 

RF n ~ 2 \ A C CP^" 2 \ A 

We obtain this submanifold by placing the IV-masses anywhere along 
the real axis 1 C C, arranged so their center of mass is zero and so 
that there are no collisions, and then taking the quotient. In other 
words, MP^ -2 \ A is the quotient of c CA -1 by dilations and 

real rehections. 


3. Main result 

In case N = 3, with the potential (eq. [Tj) above, Y 3 is a pair of pants 
- a sphere minus three points. The point of [1] was to show that the 
metric on Y 3 just described is hyperbolic provided mi = m 2 = m 3 . 
Specifically, in this equal mass case the Gaussian curvature of the met¬ 
ric on the surface Y 3 is negative everywhere except at two points (these 
being the “Lagrange points” corresponding to equilateral triangles.) 
What about Y 4 ? 

Theorem 1. Consider the Jacobi-Maupertuis metric on Y 4 induced 
as above for the case of 4 equal masses under the strong force 1/r 2 
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potential (eq. [7]). Then there are two-planes a tangent to I 4 at which 
the Riemannian sectional curvature 1C (a) is positive. 

Remark 1. The two-planes cr of the theorem pass through spe¬ 
cial points p G MP 2 C CP 2 which represent certain special collinear 
configurations. See figure [l] The two-plane a at p will be the or¬ 
thogonal complement to T p MP 2 , the normal 2-plane, and is realized as 
a = fTpMP 2 , using the standard complex structure on CP 2 . 

Remark 2. [Negative curvatures] The MP 2 of the previous remark 
is a totally geodesic surface fixed by an isometric involution. There are 
other such totally geodesic surfaces defined as fixed loci of symmetries, 
and computer experiments suggest that these all have negative Gauss¬ 
ian curvature everywhere while their normal 2-planes can have positive 
sectional curvature at some points, like our special case MP 2 . Computer 
experiments also indicate that in the direction of the normal plane there 
is positive sectional curvature over all collinear configurations of MP 2 
and not just the special configurations verified in the theorem. An 
analytic proof of these claims beyond our special case however looks 
frightening. 

Open Question. A geodesic flow can still be hyperbolic as a flow, 
without the underlying metric having all sectional curvatures negative. 
Is geodesic flow on Y 4 hyperbolic as a flow? Is it even partially hyper¬ 
bolic? 


4. Proof of the theorem 

We take the case iV = 4 in the above considerations. When all 
the masses are equal to 1 then the mass metric, used to compute the 
kinetic energy and moment of inertia, is the standard Hermitian metric 
in coordinate (<?i, Qg, <? 3 , < 74 ) G C 4 , where the g* represent the postions 
of the ith body. We reduce by translations by going to the center- 
of-mass-zero space which is a 3-dimensional subspace C 3 C C 4 having 
Jacobi coordinates as Hermitian orthonormal coordinates : 


C 3 4 C 4 given by matrix 



V2 

0 

0 


0 ' 
0 

'V2J 


in standard bases. 


As is well-known, if we start tangent to the center-of-mass-zero sub¬ 
space L( C 3 ) we stay tangent to it. Hence we can restrict the dynamics, 
potential, metric, etc. to the center-of-mass zero subspace. We de¬ 
note the potential restricted to the center of mass zero subspace in 
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Jacobi coordinates as Ul = U o L and still write ds 2 JM = U^ds 2 for the 
restricted JM metric on C 3 \A where ds 2 is the standard metric on C 3 . 

Continuing along the outline above, we now quotient by scaling and 
rotation isometries, C*, of ds 2 JM to obtain the “shape space” Y 4 and 
we label the quotient map n : C 3 \A —» I4 which takes a configuration 
q to it’s orbit C *q. We denote the vertical and horizontal distributions 

dn 

as V p = kerd p 7r = C p and % p = = T^Y^. Requiring dix\n, ds 2 JM \ n 

to be an isometry defines our induced metric on I4 whose geodesics 
correspond to iV-body motions in “shape space”. Under this induced 
metric on Y 4 we denote sectional curvature through the plane a G 
T^)Y a by /C(cr). 

Suppressing the notation of evaluating at a representative p G 7r (q), 
our main tool in the computation of JC(a), the ds 2 JM curvature, is the 
equation: 

(2) 

ultcw) = j((a 1 f/ 1 ) 2 +(a 2 c/ i ,) 2 )-||vt//2|| 2 -c/ 1 /2(a 1 2 c/ t +9 2 2 c/ L )+3Aj(« 1 

Here d a f denotes d.f(v a ) where / G C°°(C 3 ) and where a = 1,2 
with V\ , V '2 G T~L being ds 2 -or1110normal vectors whose pushforwards 
dnv a span a. The || ||, V refer to the norm, metric, and Levi-Civita 
connection for the Euclidean metric ds 2 . For the derivation of (eq. [2]) 
see Appendix A. 

The collinear configurations form a totally geodesic MP 2 C CP 2 
which is the image under the projection n of the real 2-sphere in C 3 
whose points we parameterize as 

p = (cos (j) cos 6 , cos 0 sin 9, sin 0). 

We evaluate (eq. [2]) and find positive sectional curvature over the 
configurations with 9 = 7t/2 (see figure [lj in the direction of the zTIRP 2 
plane. This plane is spanned by the pushforwards of 

dp 

Vi = —i— = i (sin 0cos0, sin 0 sin 6>, — cos0) 

<70 

v 2 = 1 = i(— sin 9, cos 6, 0). 

cos 0 o0 

TERMS 1: Over MP 2 in the iTMP 2 direction, the last term and 
first two terms (the first partials) of (eq. [^) vanish: 

V] ■ iv 2 = 0, d a U L = 0. 

That V\ ■ iv 2 = 0 is clear: i rotates v 2 into purely real coordinates. To 
evaluate the 1st partials, note Lp has purely real coordinates and VU 
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2 4 3 1 

H— • I X I • II * 

Center of Mass 


FIGURE 1. : The collincar configurations p which we consider, 
has kth component c/j d <?fc ; so V| l p U has purely real coordinates. 

' jk 

Now since Lv a has purely complex coordinates, 
d a U L = V\ Lp U ■ Lv a = 0. □ 


TERMS 2: With the notation Lp = (q 1 , q 2 , q$, qf), Lv a = i(u 4 , v 2 , v%, u 4 ) 
and pj k = q .\ qk , ctjk = {y{ — v\) 2 + (v 2 — v 2 ) 2 G M, the 2nd partials 
terms of (eq. [![) are: 


dlU L + d 2 2 U L = -2Y j a ]k p%. 

j>k 

We write our standard coordinates on C 4 as qj = x 3 +iyj, then since 
Lv a is purely imaginary: 

dU d 2 U 

alu L = vi^vtr ■ lv .) • L Va = (v| ir ^) • Lv. = ^ 

d 2 U d 2 U 

Next we compute Q \ Lp = 2 p 4 fc for j ^ k and ~q^\lp = ~ 2 J2p% 


\L P v k a vl 


dyjdVk 


3+k 


so now: 


dpL = 


- 2 Y,P%( C 

j¥=k 


,k\2 


■v’A) 


-2 E4(M) 2 -2»„v.+(»;) ! ) = -2£ 

j>k j>k 


□ 


RESULT: Over t/ie circ/e 6 = n/2, /C(iTIRP 2 ) is positive. 

For, substituting terms 1 and 2 into formula (eq. [2]), we see that: 


0 < UlK. = —1| VU/2|| : 


+ Ul^ 2 a jk p% 

j>k 


0 < K, 
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(3) E(E p %) * 2 < (E p%) (E a Jkp%) 

k j^k j>k j>k 


Taking 9 = 7 t/2 and with the notation introduced in terms 2, we find 
the relations: 


Pn v^COS 4>' Vasili. 


Pl3 — 


y/2 


p 14 — 


cos <p — Sill < 

V2 

cos 0 + sin < 


— — P 24 


— —p 23 


1 1 

12 — r; a 34 — —; t~ 
P34 Pi 2 

1 

T13 — o-f 1 = a 24 

Pl4 

1 

Cl 14 = “5-h 1 = Cl23- 

PlS 


Now the left side of (eq. [3]) works out to: 

2 ((p ?2 + P 13 + Pu ) 2 + (Pi 3 — P 14 — P 34 ) 2 ) = 


2 S Pifc+ 2 Pl2(Pl3+Pl4)+ 2 P34(Pl4 P 13 )) - 2 E P% 96 gj n 2 20 COS 2 20 


k>j 


k>j 


= 2 ^ p® k + negative term 


k>j 


and the right side of (eq. [3]) works out to: 

(Pl 2 + P 34 + 2 (Pl3 + P 14 )) ( + 2 (Pl3 + Pl4 + + tt)) = 

P 34 P 12 Pl4 Pl3 


8 -)(sin 2 20(p^ 2 +P3 4 ) + CQS o 2 < ^ (Pi 3 +Pi 4 )+ 2 (Pi 3 +Pi 4 ) = 


"sin 20 cos 2 20 " 


4 16 

2 E P %+ cot2 2 < ^(Pi 3 +Pi 4 )+8 tan 2 2 0 (Pi 2 +P 3 4 ) + (P? 3 +Pi 4 )( 2 9 , + 2 9 a 

sm zcp cos zcp 
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= 2^2 pj k + positive term. 

k>j 

Therefore the inequality (eq. g holds! □ 
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Appendix A. DERIVATION OF eq. [2} 

Take a ds 2 -orthonormal basis {u a } for C 3 with V\,v 2 G R p . 

The Kulkarni-Nomizu [K-N] product formula for conformal curva¬ 
tures (0, pg. 51) reads: 

Rabcd U L R abcd \^ds j dlL dll ® dll ~\~ -\\du\\ ds ^)}abcd 

where u := | log Ul and the overbars denote curvature with respect to 
the ds 2 JM - metric and all other quantities (no overbars) are with respect 
to the fi.s' 2 -metric. Then R ab cd = 0 since ds 2 is the flat Euclidean metric 
of C 3 = M 6 . Taking cd = ab we have: 

U\K ab = R abab = -U L (’Vdu bb +Vdu aa -du b ®du b -du a ®du a + \\du\\ 2 ) = 

= -U L (d 2 a ii + d 2 b u - ( d a uf - (d b u) 2 + || Vu|| 2 ). 

Next O’Neill’s formula ([2], pg. 213) gives 


/C(d7TUi, dTTV 2 ) = K 12 + U 2] 1 \ 2 ds 2 jM 


where V a = 

We then compute: 



and X v denotes ds 2 JM projection of X onto V. 


_ d a U L _ V | Lp U ■ Lv a 
Z 2U L ~ 2U L (p) 


2 _ dlU L {d a U L f _ V|Lp(W • LVg) • LVg (dgU L ) 2 
aU 2 U L 2 U 2 2 U L (p) 2 U L ( P y 


and 
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Note that XU G {q G C 4 : Qj — 0} and Lv a is a ds 2 orthonormal 
basis for this center of mass zero subspace, hence 

l|Vf/|| 2 = £(W . Lp„) 2 = Y.(W 2 = iulw V„|| 2 . 

Substitution into the K-N formula gives 


(4) k 12 = -^Y (a l UL + ^ ~ i (SlC/ " + d2U ^ + ll vc/ / 2 ll 2 )- 

To compute O’Neill’s Lie bracket term we write our standard coor¬ 
dinates on C 3 as ( x 1 + ix 2 , ...,x 5 + ix 6 ). 

Let H i = X 3 d x j , H 2 = Y 3 d x j G ^ be any horizontal vector fields. 
The vertical vector fields are spanned by the Euler vector field E = 
x 3 d x :i and iE. Then Hj ■ E — Hj ■ iE = 0 and: 

[//,, H 2 )-E = Y,X j x k d xj Y k - Y 3 x k d x: ,X k = 

k 




= £ X j (x k Y k )-5 k Y k )-Y j (d xj (x k X k )-5 k X k ) = £ 

fc k 

and likewise: 

[Ri,R 2 ]-hE = £ {Y 3 d x jX k —X 3 d x jY k )x k+1 +(X 3 d x jY k+1 —Y 3 d xJ X k+1 )t 

fc odd 

= 2 ^ + X k+1 Y k = 2 Hi ■ %H 2 . 


k odd 


Then 


IK,Upl 2 = ds) M (K,vy, , ) 2 +<fe 2 «(K,VI, 


iEr, 


\pWUl(p)' 


\pWu l (p)‘ 


ul 


-m,v 2 ] -Emv h v 2 h E ?) = LLLdLilili - 


IpI^l" 1 ’ 1 ’ ’ 2J “ 2 ' u "’ ' 21 ' |p| 2 I4(p)br 

Now substitution of this Lie bracket expression and (eq. [f]) into 


O’Neill’s 


uuotituiiuii uiiio i_jic uiciuivcl CApiraDiuii ciiiu. 

formula and multiplying by U\ yields (eq. [2|. □ 
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